Renewal Aging as Emerging Property of Phase Synchronization 
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In this letter we examine a model recently proposed to produce phase synchronization [K. Wood et 
al, Phys. Rev. Lett. 96, 145701 (2006)] and we show that the onset to synchronization corresponds 
to the emergence of an intermittent process that is non-Poisson and renewal at the same time. We 
argue that this makes the model appropriate for the physics of blinking quantum dots, and the 
dynamics of human brain as well. 
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Phase synchronization of coupled clocks (oscillators) is 
a growing field of research, which is fast developing from 
the seminal work of Winfree 0] and Kuramoto 0] . Using 
the model of coupled clocks of these authors it is possible 
to derive yj the Turing structure which, in turn, trig- 
gered the research work in the field of diffusion reaction 
p. Physiologists are using a clock as a representation 
of a single neuron. In the work of Ref. ^ a neuron is a 
clock whose behavior is described by a chaotic attractor: 
the Rossler oscillator . The authors of Ref. (Sj and of 
Ref. (21] have shown that coupled stochastic clocks can 
show cooperative (synchronized) behavior. Brain func- 
tions, such as a cognitive act, rest on the cooperative 
behavior of a collection of many neurons The au- 

thors of Ref. study the dynamic of the cooperation 
of a collection of many neurons by mapping the brain ac- 
tivity into a network. They find that the changes in the 
topology of the network describing the brain activity are 
driven by a non-Poisson renewal process operating in the 
non-ergodic regime 12] . A collection of Blinking Quan- 
tum Dots (BQDs) |1^ has the same dynamical property 
as the brain activity. The non-Poisson non-ergodic char- 
acter of the distribution of the sojourn times of the BQDs 
in the "light" and in the "dark" state is a well estab- 
lished property . The renewal character of the BQDs 
dynamic has been established only recently [l^ . 

The similarity between the BQDs and the brain activ- 
ity dynamics makes it plausible to search for a dynamic 
model accounting for both complex systems. In this Let- 
ter we show that a simplified version of the model of 
ReL Q affords significant suggestions on how to realize 
this important purpose. The authors of ReL use a sys- 
tem of coupled 3-state stochastic clocks, while we use a 
system coupled 2-state stochastic clocks. We shall prove 
that at the onset of phase synchronization the dynam- 
ics of the system has the same properties as the BQDs 
[T3L 113 ] and brain dynamics ■ They are non-Poisson 
renewal processes operating in the non-ergodic regime. 

We consider a Gibbs ensemble of systems with N 
2-state stochastic clocks, each of them coupled to Nc 
clocks. We denote by ]1> and |2> the two states of 
the clock, corresponding to the phases ^—0 and (f> 



respectively. The master equation for a single clock of a 
system of the Gibbs ensemble is 
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Pi {P2) is the probability of finding the clock in the state 
]1> (12>), and 1712 (<?2i) is the rate of transitions from 
the state ll> (12>) to the state |2> (]!>). The transition 
rates gi-j^and §21 are defined by means of the prescription 
of Ref. 



512(21) = 9 exp [K{-K2(1) - 71-1(2))] 



(2) 



In Eq. , .9 is the unperturbed transition rate of a single 
clock, A'>0 is the coupling constant and tti (7r2) is the 
fraction of the Nc coupled clocks that are in the state 1 1 > 
(|2>). The authors of ReL place their clocks in a d- 
dimensional lattice where only the nearest neighbors are 
coupled. Thus, every clock is coupled to Nc=2d clocks. 
We adopt, instead, an all to all couphng: Nc=N—l. 
With this choice, when N—^oo, the mean field approx- 
imation I15I 
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is an exact property. Thanks to the mean field approx- 
imation of Eq. Q and to the normalization condition 
P1+P2—I, the master equation Eq. reduces to 



dU{t) 
dt 



-2g cosh( j^:n)n + 2g sinh(i^:n) = - ^ gj^'' , (4) 



where n=Pi-P2 and ne[-l,l]. Eq. describes the 
overdampcd motion of a particle, whose position is IT, 
within the potential ^(n) [l|]. Using Eq. (g)), we find 
that the potential 1^(11) is symmetric and the values of its 
minima depend only on the coupling constant K. More- 
over, we find that there is a critical value Kc of the cou- 
pling parameter K 



K = Kr 



(5) 



such that: 1) If K<Kc the potential V^(n) has only one 
TT, minimum for n,iiin=0; 2) if K>Kc the potential V^(n) 
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is symmetric and has two minima ibllmin separated by a 
barrier with the maximum centered in 11=0. As shown in 
Fig. n the value n,„in and the height of the barrier VfO ) 
are increasing function of the coupHng constant K |l7| . 
In particular nmin^l and ^(0)— >+oo when if — >+oo. 
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FIG. 1: The potential V'(n), rescaled by a factor 10 , as a 
function of 11 for ^=0.01 and different values of the coupling 
constant K. _K'=1.05 for the continuous line, 7^=1.1 for the 
dashed line, and K=1.2 for the dot-dashed line. 



The time evolution of the variable 11 is determined by 
the minima and maxima of the potential V(n). Thus, 
two kinds of dynamical evolution are possible: 1) If 
K<Kc, n(t) will reach, after a transient, an asymptotic 
value n(cx3)=0 not depending on the initial conditions 
n(0); 2) if K>Kc, n(t) will reach, after a transient, an 
asymptotic value U{oo)=+Ilnun^O (n(oo)=— IIininT^O) 
for an initial condition n(0)>0 (n(0)<0), while the ini- 
tial condition n(0)=0 wiU resuh in n(t)=0 Vi. In Fig. El 
we compare the minima ibllfnin of the potential T^(n) and 
the numerical evaluation of n(oo) for different values of 
the coupling constant K. Fig. [21 shows that a phase tran- 
sition occurs at K=Kc=l- For a single clock the 
condition n(oo)=±ni„in=-Pi— -P27^0 corresponds to the 
statistical "preference" to be either in the state |1> or 
|2>. This is a consequence of the fact that the transi- 
tion rates gi2 and g2i of Eq. |2Il are different if K>Kc- 
Plugging Eq. ^ into Eq. Q and allowing 11 to reach its 
asymptotic value, we get 

gi2 ^ gexp{-KU{oo)) ^ g2i ^ g exp{KIl{oo)) . (6) 

Fig. O confirms the prediction of Eq. © showing that 
if n(oo)=-|-n„iin (II(cx3)=— Ilinin) the siugle clock spends 
on average more time in the state |1> (|2>). The proba- 
bility density function for the sojourn times in both the 
preferred and not-preferred state are exponential func- 
tions with different mean sojourn times. 

Let us now explore the collective behavior of a single 
system of N clocks under the all to all coupling condition. 
Following the authors of Ref. 9], we define the global 
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FIG. 2: The minima illmin of the potential V^(n), and the 
asymptotic value n(oo) as a function of the coupling constant 
K. The full line is the theoretical prediction obtained solving 

^ =0 and ^-^^>0. The squares denote the result of 
the numerical evaluation of n(cx3) with a Gibbs ensemble of 
systems each with Ai'=10000 clocks. 
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FIG. 3: The typical time evolution of a single clock of a sys- 
tem in the Gibbs ensemble. The solid (dashed) line refers 
to the case n(oo)=-|-nmin (n(cx3)=— flmin). The number of 
clocks of the single system of the Gibbs ensemble is 10000, 
the unperturbed transition rate is g—0.01 and the coupling 
constant is K—1.05. 
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The symbol i indicates the imaginary unit, is the 
phase of the j-th clock, (tt) if the clock is in the state 
|1> (|2>), and Ni{t) (iV2(i)) is the number of clocks of 
the system in state |1> (|2>) at a time t. When iV— >oo, 
the single system becomes a Gibbs ensemble on its own 
as all the clocks are identical, and, at the same time, the 
mean field approximation (Eq. (Q) becomes valid. In 
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this case the master equation for the single system is the 
master equation Eq. where Pi {P2) is now the prob- 
abihty of finding a clock of the system in the state |1> 
(|2>). Finally, from Eq. {Tj), in the limiting case N^oo, 
we get 



(8) 



Thus, for a system with infinitely many clocks, the time 
behavior of the global clock variable ^ is identical to that 
of the variable 11 of Eq. for K>Kc (we exclude the 
initial unstable condition ^(0)=0) 



^(00) = ±n„,in ^0 if m k 0. 



(9) 



The condition ^(00)7^0 of Eq. ^ proves that a global 
phase synchronization occurs in the system at the onset 
of phase transition {K>Kc). Moreover, the time evolu- 
tion of a single clock of the system is the one depicted 
by Fig. O state |1> (|2>) is statistically preferred if 

f (00) =+ninm {^{oo) = -Umin). 

Is the phase synchronization of Eq. lO present in a sys- 
tem with a finite number of clocks? From the definition 
of transition rate and of probability, it follows that 
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where Ni^2 {^2^1) is the number of clocks that undergo 
a transition from state |1> (|2>) to state |2> (|1>) per 
unit of time, and N is the number of clocks of the system. 
Using the law of large numbers we get that, for any 
N fitting the condition oo>iV::^l, 



-j^ = 912P1 

— .921^2 



£12^1 

£21^2 



(11) 



where £12 and £21 are fluctuating variables whose intensi- 
ties are ccl / y^N) . From Eqs. (|10(l and we conclude 
that the master equation of a system with a finite number 
of clocks is equivalent to that of a system with an infinite 
number of clocks whose transition rates fluctuate: 




-(.gi2 + £i2)Pi + (.921 + £21)^2 
-(.921 + £21)^2 + (.912 + £l2)Pl 



(12) 



If oo>A''^l, we can still consider the mean field approxi- 
mation of Eq. to be valid. Using the master equation 
Eq. (|12|) and the normalization condition Pi+P2 = l, we 
get for the variable n=Pi— P2 the following equation of 
motion 



dn{t) dV{Il) 



dt 



an 



7im{t) + B{t). 



(13) 



The presence of the fluctuations ?7=ei2 + £21 and 9~ei2 — 
£21 in Eq. H13|l has the effect of triggering transitions from 
one well to the other of the potential V{Ii) of Fig. ^ 



Thus, for a system with a finite number of clocks the 
phase synchronization of Eq. is not stable. The global 
clock variable ^ of Eq.(f7|) fluctuates (for K>Kc) between 
the two minima, ibllniin, of the potential 1^(11), as con- 
firmed by Fig. 2] The single clock follows the fluctuations 
of the global clock variable ^, switching back and forth 
from the condition where the state 1 1> is statistically pre- 
ferred (time evolution described by the full line of Fig.O 
to that where the state |2> is (time evolution described 
by the dashed fine of Fig. 13). 
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FIG. 4: The global variable ^{t) as a function of time, for 
K=1.05, g—0.01 and a system of 1000 clocks. 

The probability density functions of the sojourn times 
in the state ^>0 or ^<0 (Fig. ^ are identical since the 
potential 1^(11) of Fig. ^ is symmetric. Thus, we de- 
note them with the same symbol iP{t). Let us consider 
a condition where the coupling constant K is close to 
the critical value K—Kc'^1- In this case, the height of 
the barrier V{0) dividing the two wells of the potential 
y(n) of Fig.n]is smaller than or comparable to the inten- 
sity of the fluctuations?? and 9 of Eq. . Under this 
condition, we expect |l3| 'il^{T)(y:l/T^-^ for an extended 
interval of sojourn times. This is exactly what we ob- 
serve in Fig. [SI where the full line denotes the survival 
probability ^'(t), namely, the probability of observing 
a sojourn time larger than r pO]. For any fixed value 
of the unperturbed rate g (Eq. jSJl), the height V{0) of 
the barrier dividing the two wells of the potential 1^(11) 
(Fig. ^1 increases as the value of the coupling parame- 
ter increases. Eventually, the height V{Q) will become 
much larger than the intensity of the fluctuations ry and 
6 (Eq. (|13|l l. As a consequence, the power law behavior 
('0(T)ocl/r^-^) of Fig. 13 disappears, the theoretical argu- 
ments of [T9I I loses validity, and an exponential behavior 
emerges, as predicted by the Kramers theory jld] . 

Finally, we show that the transition between the states 
^>0 and ^<0 shown in Fig. ^ is a renewal process. For 
this purpose we use the aging experiment of Ref. |l4j . 
We evaluate the survival probability 'i'{ta,T) of age ta- 
This is the probability of observing a sojourn time larger 
than T if the observation starts at a time ta after a cross- 
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FIG. 5: The survival probability 'l'(r) (full line), the survival 
probability *l/(ta, r) (dashed line) and the survival probability 
in the renewal case ^r(ta, t) (full thick line) as a function of 
sojourn time r. Here, we plot for clarity only the survival 
probability of age ta=500. 

ing from ^>0 to i^<0 or vice versa {'^/{0,t)—'^{t)). We, 
then, compare '^(ta,T) with the expected survival prob- 
abihty 'i'r(ta,'T) of age ta, evaluated according to the re- 
newal theory If ^(i^, T)=^'r(ta, ''') the process 
described by the survival probability v1/(t) is a renewal 
process. If \E'(r) is not an exponential function, 'i/{ta,T) 
yields a slower decay than ^'(r), a condition that is de- 
noted as "aging" jlj| . Fig. El shows that the transition 
between the states ^>0 and ^<0 (Fig. ^ is a renewal 
process and that there is aging. With increasing values 



of coupling parameter K the renewal property is not lost, 
but the aging property is because the survival probability 
^'(r) becomes an exponential function 

In conclusion, we have shown that, for a Gibbs en- 
semble of systems with N=oo coupled 2-state clocks, a 
phase transition occurs at the critical value of coupling 
parameter Kc—1 (Fig. I^J- The phase transition mirrors 
a statistical "preference" for a clock of one of the Gibbs 
ensemble systems to be in one of the two possible states 
(Eq.®). For a single system with N=oo coupled 2-state 
clocks, the phase transition of Fig. [5] signals the onset of 
a global phase synchronization (Eqs. O, © and @). 
If the number of clocks of the system is finite the phase 
synchronization is not stable. In this case, the variable 
^ (Eq. [7| describing the collective motion of the system 
(Fig. ^ is characterized by non-Poisson intermittent be- 
havior. At the onset of the phase transition (K—Kc-^l), 
the zero crossings of the global variable ^ (Fig. ^ define 
a series of events with the following properties: 1) The 
probability density function of inter events intervals (full 
line of Fig. O has a non-Poisson non-ergodic character 
{'iJj{t)(x1/t^'^^ infinite mean inter event interval). 2) 
The sequence of inter events intervals satisfy the renewal 
aging condition (Fig.EJ. The properties 1) and 2) are ob- 
served properties of the events in both BQDs [isL IT3. IToj l 
and brain activity , suggesting that a system consist- 
ing of a finite number of coupled 2-state clocks may be a 
good model for the dynamics of both processes. 
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